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Motivation 

•  Urgent  need  for  effective  padding  materials  to 
minimize  severity  of  blast-related  injuries  due  to 
improvised  explosive  devices  (lEDs) 

•  For  use  in  helmets,  personnel  protective  equipment 
(PPE),  vehicle  interiors,  etc. 

•  Difficulty  in  injury  detection  and  the  resulting 
intractability  of  injuries  results  in  long-term 
treatment,  suffering,  and  costs  to  the  Army 

•  Need  a  method  for  rapidly  evaluating  efficacy  of 
candidate  materials  for  various  applications 
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Determination  of  Viscoelastic  Properties  of 
Foam  Materials 

•  GOAL:  Determine  material  response  (transfer  function)  in  the  frequency  domain 

•  RATIONALE:  Preliminary  part  of  blast  wave  imparts  an  impulsive  force  on  the 
victim 

•  Input  provided  by  impact  hammer;  force  =  F(s) 

•  Output  is  the  acceleration  of  the  mass,  A(s) 


•  Fractional-order  models  that  feature  a  three  parameter  model  in  which 
"fractional  derivatives"  PROVIDE  MORE  ACCURACY  AND  FLEXIBILITY  IN 
MODELING  DYNAMIC  BEHAVIOR  OF  FOAM. 
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Challenges 

•  Most  padding  materials  exhibit  viscoelastic  behavior 

•  Dynamic  behavior  is  nonlinear  and  frequency 
dependent 

•  Cannot  be  accurately  modeled  using  simple  lumped- 
parameter  elements 

-  spring,  k  (energy  storage) 

-  damper,  b  (energy  dissipation 
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Determination  of  Viscoelastic  Properties  of 
Foam  Materials 

•  Previously  considered  two  approaches  to  modeling  dynamic  behavior  of  foam 

1.  Integer-order  transfer  function,  featuring  "traditional"  linear  springs  and  dampers, 
i.e.,  Hooke's  Law  springs  and  viscous  damping,  respectively 

2.  Fractional-order  models  that  feature  a  three  parameter  model  in  which  "fractional 
derivatives"  are  used 


•  SIMULATIONS  INDICATE  THAT  OPTION  #2  PROVIDES  MORE  ACCURACY  AND 
FLEXIBILITY  IN  MODELING  DYNAMIC  BEHAVIOR  OF  FOAM. 
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Fractional  Derivative??? 


a  =  0 


a  =  1 


spring 

F  =  kx 


viscoelastic  material 

F  =  yx(a) 


damper 

F  =  bx 
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Fractional  Derivatives  to  Model  Foam  Behavior 


•  Use  a  three-parameter  impedance-based  (frequency  domain)  model  to 
represent  viscoelastic  behavior  of  the  foam  based  on  force  input  and 
acceleration  response  data 

■  G0  =  shear  modulus 

■  Gx  =  modulus  that  represents  shear  stiffness  and  dissipative  behavior  of  the  foam 

■  a  =  order  of  the  fractional  derivative;  used  to  determine  G0  and  Gx 


represents  pure 
stiffness  behavior  of 
material 


represents  a  combination  of 
stiffness  and  dissipation 
behavior  of  material 


G0+saG\  = 


model 


2  17^1^  1 

V  F(s)J 

experimental  data 


mt 
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Advantages 

•  Ability  to  use  lumped-parameter  modeling 
techniques  for  simulations 

-  simplified;  increased  accessibility 

-  more  unified  modeling  approach;  fewer  interfaces 
between  models 
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Overall  Modeling  Approach 


CFD  model  provides 
incident  blast  pressure 
time  history 


bond  graph  model  with 
integer-order  and 
fractional-order  lumped 
parameter  elements 
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Padding  Materials  Used  in  Testing 


•  Used  two  types  of  commercially  available  foam  in  lieu  of 
candidate  padding  materials  to  be  provided  by  ARL 


padding  material 

density 

(g/cm3) 

cell 

type 

color 

Ethylene  Vinyl  Acetate 

0.032 

closed 

white 

Acrylonitrile  Butadiene  Rubber  / 
Polyvinyl  Chloride 

0.088 

closed 

black 
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Determining  Model  from  Experimental  Data 


1.  Acquire  input  force  and 
acceleration  response  data  from 
impact  testing 

2.  Compute  averaged  power 
spectra  from  five  experimental 
runs 

3.  Apply  curve-fitting  algorithm  to 
compute  a  plot  of  the  transfer 
function  between  input  force 
and  acceleration  of  mass 
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Determining  Model  from  Experimental  Data 


4.  Compute  real  and  imaginary 
parts  of  transfer  function  data 


frequency  (Hz) 


5.  Perform  a  linear  regression  on 
each  part 
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Transfer  Function  Curve  Fit  Results 


cm  40 

CO 

1 

Z 

experimental  data 
constant  GQ 

f— 

L 

aT  20 

"n 

_ G0=G0(m) 

u 

3 

■M 
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'c 

CD 

E  u0 

100  200 

300 
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Results:  Material  Parameters 


padding  material 

color 

Q 

Go  (N/m2) 

Gi  (N-sQ/m2) 

Ethylene  Vinyl  Acetate 

white 

0.56 

3.00  x  10s 

3.19  x  103 

Acrylonitrile  Butadiene  Rubber  /  PVC 

black 

0.60 

2.39  x  10s 

3.05  x  103 
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Validation  of  Model  vs.  Experimental  Data 

Ethylene  Vinyl  Acetate  (white  foam) 
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Material  Model  Simulations 


GOALS 


1.  Verify  that  three-parameter  model  material  properties  that  were  determined 
from  impact  tests  are  accurate 

2.  Determine  onset  of  saturation:  material  can  no  longer  store  or  dissipate  energy 
due  to  blast  pressures;  maximally  compressed  during  loading;  acts  like  a  mass, 
i.e.,  compliance  and  dissipative  properties  are  effectively  equal  to  zero 


helmet  shell  material 
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Blast  Pressure  Time  Histories  from 
CFD  Model 
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Model  Assumptions 

•  Frontal  area  of  helmet  =  0.05  m2 

•  Foam  pad  thickness  =  5mm 

•  Saturation  occurs  when  pad  is  compressed  down  to 
1mm 

•  Victim  is  knocked  over  by  blast,  i.e.,  no  impact  with 
objects  such  as  a  vehicle,  structure,  etc. 

-  results  in  a  free-free  boundary  condition  for  helmet  mass 
and  mass  of  the  victim's  head 

-  victim's  resistance  to  motion  can  be  easily  incorporated 
into  model 
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Example  Simulation  -  Saturation 
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Summary  of  Dynamic  Modeling  Results 


foam 

blast  ID 

peak  blast 
pressure 

peak  blast 
pressure 

peak  blast 
force 

peak  blast 
force 

peak 

transmitted 

force 

peak 

transmitted 

force 

attenuation 
of  peak  force 
(trans/blast) 

max 

compression 

saturation 

Pa 

lbf/inA2 

N 

Ibf 

N 

Ibf 

-- 

mm 

white 

1 

1.57E+04 

2.27E+00 

9.91 

2.23 

4.53 

1.02 

46% 

0.013 

NO 

black 

1 

1.57E+04 

2.27E+00 

9.91 

2.23 

4.22 

0.95 

43% 

0.014 

NO 

white 

3 

3.28E+04 

4.76E+00 

21.1 

4.74 

9.77 

2.20 

46% 

0.028 

NO 

black 

3 

3.28E+04 

4.76E+00 

21.1 

4.74 

9.16 

2.06 

43% 

0.028 

NO 

white 

7 

8.60E+04 

1.25E+01 

55.5 

12.50 

24.6 

5.53 

44% 

0.071 

NO 

black 

7 

8.60E+04 

1.25E+01 

55.5 

12.50 

23.1 

5.19 

42% 

0.071 

NO 

white 

25 

3.15E+05 

4.57E+01 

203 

45.70 

90.6 

20.37 

45% 

0.261 

NO 

black 

25 

3.15E+05 

4.57E+01 

203 

45.70 

85 

19.11 

42% 

0.267 

NO 

white 

31 

4.10E+05 

5.95E+01 

265 

59.60 

117 

26.30 

44% 

0.339 

NO 

black 

31 

4.10E+05 

5.95E+01 

265 

59.60 

110 

24.73 

42% 

0.341 

NO 

white 

125 

3.25E+06 

4.71E+02 

2037 

457.00 

376 

84.53 

18%] 

1.000 

YES 

black 

125 

3.25E+06 

4.71E+02 

2037 

457.00 

368 

82.73 

1.000 

YES 
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Summary 


•  Fractional-order  modeling  of  material  models 
provides  accurate  results 

•  Modeling  technique  provides  access  to  multiple 
metrics  that  can  be  used  to  quantify  desired  padding 
material  properties 

•  Overall  approach  can  be  expanded  to  include  wider 
range  of  blast  scenarios  and  PPE  applications 
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Future  Work 


•  Improve  computation  performance  of  fractional 
derivative  algorithm 

•  Use  actual  helmet  and  padding  materials 

•  Expand  blast  scenarios  beyond  helmet  applications 

-  vehicle  interior  padding 

-  other  PPE 

•  Investigate  other  dynamically  complex  materials  that 
are  of  interest  to  the  Army 

-  crushable  materials 

-  others??? 
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This  paper  concerns  the  development  of  a  combined  experimental  and  simulation  approach  to  assist  in  rapidly 
quantifying  metrics  that  characterize  the  energy  absorption/ dissipation  capacity  of  candidate  materials  to  be 
used  in  retrofitted  helmets  and  other  personnel  protective  equipment  (PPE)  in  an  explosion  environment. 
The  dynamics  of  a  cross  section  of  a  helmet,  padding  material,  and  a  proportional  mass  that  represents  the 
head  of  the  victim,  are  modeled  using  impedance-based  bond  graphs  that  feature  primitives  whose  constitutive 
behaviors  are  described  using  fractional  derivatives.  Force  inputs  to  the  model  are  derived  from  computational 
fluid  dynamic  (CFD)  models  of  improvised  explosive  devices  (IED)  with  a  range  of  net  explosive  weights 
(NEW.)  Metrics  such  as  the  saturation  point  of  the  material  with  respect  to  maximum  blast  loads,  and  the 
ratio  of  the  total  blast  energy  to  the  energy  that  is  transmitted  to  the  victim,  are  used  to  rate  two  candidate 
materials. 


1  Introduction  and  Thesis 

The  use  of  energy- absorbing/dissipating  materials  in  personnel  helmets  to  reduce  the  effects  of  blast  pressures 
can  have  unanticipated  adverse  effects.  The  presence  and/or  configuration  of  these  materials  can  focus 
energy  in  such  a  way  that  it  can  cause  unforeseen  traumatic  brain  injuries  (TBI).  There  is  a  wide  range  of 
candidate  materials,  particularly  viscoelastic  materials,  whose  abilities  to  both  dissipate  and  absorb  energy 
vary  considerably.  These  abilities  are  difficult  to  quantify  and  model  using  linear  lumped-parameter  modeling 
techniques  and  bulk  material  properties.  This  challenge  is  further  complicated,  given  the  extreme  nature  of 
the  input  forcing  functions  encountered  in  a  blast  environment,  i.e.,  short-duration  impulses  that  result  in 
high  rates  of  energy  transfer  to  a  victim. 

The  authors  posit  that  using  fractional-derviative-based  material  models  for  viscoelastic  padding  ma¬ 
terials,  allows  the  use  of  a  lumped-parameter  modeling  approach,  which  will  yield  more  accurate  results 
than  those  with  linear  material  models.  This  approach  eliminates  costly  and  time-consuming  computations 
associated  with  finite  element  modeling  approaches,  which  are  not  necessary  for  significant  portions  of  the 
model  that  are  relatively  incompressible,  i.e.,  the  helmet  shell  material  and  the  victim’s  head,  as  compared 
to  the  helmet  padding  material. 

2  Overall  Approach 

The  overall  approach  of  the  work  presented  in  this  paper  is  shown  in  Figure  1.  To  assess  the  dynamic 
behavior  of  the  helmet  padding  material  applied  as  a  protective  material  in  explosive  environments,  it  was 
necessary  to  determine  input  forces  that  are  imparted  to  a  victim  in  a  wide  range  of  explosive  events,  i.e., 
of  varying  intensity.  Clutter  and  Connolly  [2]  performed  an  interdisciplinary  study  in  which  computational 
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CFD  blast  simulations  that  include  a  "virtual  victim" 


experimentally  determined  material  parameters 
via  auto  &  cross-power  spectral  density  function  analysis 
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Figure  1:  Overall  Approach 


fluid  dynamic  analysis  (CFD)  was  used  to  generate  time  histories  of  blast  pressures  as  experienced  by  a 
virtual  victim.  These  time  histories  were  used  to  compute  Input  forces,  as  determined  by  the  pressures  at 
the  center  of  the  victim’s  helmet,  as  shown  in  Figure  2.  An  example  of  an  input  force  time  history  for  a 
particular  net  explosive  weight,  i.e.,  explosion  intensity,  is  shown  in  Figure  3.  Material  properties  of  the  foam 


Figure  2:  CFD  model  of  blast  interaction  with  victim 
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padding  materials  are  determined  by  subjecting  them  to  a  impulsive  force  input,  which  is  similar  to  that  of 
an  explosion,  albeit  on  a  smaller  scale.  Acceleration  and  force  data  from  these  experiments  are  used  to  form 
experimentally  determined  transfer  functions  -  that  feature  fractional  exponents  of  the  complex  frequency,  s, 
from  which  material  properties  are  extracted.  The  experimental  setup  and  protocols  are  detailed  in  Section 
3  and  the  procedure  for  determination  of  material  properties  in  the  context  of  a  fractional-derivative  model 
are  detailed  in  Section  4. 


3  Experimental  Setup  and  Protocols 

The  experimental  test  setup  is  shown  in  Figure  4.  The  test  article  consisted  of  two  pieces  of  viscoelastic 
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Figure  4 ■  Experimental  setup 


foam,  0.25  inch  thick  by  1  .25  inch  square,  adhered  to  an  aluminum  block,  0.5  inch  thick  by  1.25  inch 
square,  using  a  thin  layer  of  ethyl  cyanoacrylate  applied  with  a  brush.  A  piezoelectric  accelerometer  was 
attached  in  the  z  coordinate  through  a  threaded  hole  at  the  bottom  of  the  aluminum  block.  The  test  article 
was  secured  using  a  vice  grip  and  held  in  minimal  compression.  An  impulse  hammer  implemented  with  a 
piezoelectric  accelerometer  was  used  to  tap  the  specimen.  Data  was  collected  at  a  rate  of  1  kHz  with  a  USB 
data  acquisition  accessory  and  Lab  View™  running  on  a  notebook  computer.  The  experimental  procedure 
was  repeated  five  times  for  each  type  of  padding  material. 

Issues  that  were  encountered  and  their  subsequent  resolutions  included: 

1.  A  60  Hz  ground  loop  that  was  eliminated  by  creating  one  common  ground  among  all  electronics. 

2.  Double- hit  phenomena  were  encountered  when  using  harder  tips  because  of  the  deformation  of  the 
aluminum.  Therefore,  it  was  necessary  to  use  the  softest  tip  available  in  the  kit. 

3.  Aliasing  issues  were  a  problem  due  to  the  short  nature  of  the  test,  which  were  remedied  by  using  a 
higher  sampling  rate. 

4.  The  natural  frequency  of  the  system  was  too  low  to  obtain  viable  acceleration  data  -  this  required 
using  an  aluminum  block  instead  of  titanium,  which  had  been  the  preferred  metal,  due  to  its  lower  loss 
coefficient.  Further,  this  also  required  a  thinning  of  the  foam  pads  to  increase  the  effective  stiffness  of  the 
test  article.  Two  test  specimens,  each  with  a  different  viscoelastic  padding  material  were  used  for  the  tests. 
The  physical  characteristics  of  these  materials  are  given  in  Table  1. 

4  Helmet  Padding  Model  Using  Fractional  Derivatives 

A  schematic  that  represents  a  cross  section  of  the  helmet,  padding  material,  and  a  proportional  mass  that 
represents  the  victim’s  head  is  shown  in  Figure  5.  The  corresponding  bond  graph  model  of  the  dynamics  of 
this  systems  is  shown  in  Figure  6. 
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padding  material 

density 

(g/cm3) 

cell 

type 

color 

Ethylene  Vinyl  Acetate 

0.032 

closed 

white 

Acrylonitrile  Butadiene  Rubber  / 
Polyvinyl  Chloride 

0.088 

closed 

black 

Table  1:  Padding  Material  Physical  Characteristics 
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Figure  5:  Helmet  cross  section 
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Figure  6:  Impedance-based  bond  graph  with  viscoelastic  primitive  element 


The  experimentally  determined  transfer  function  is  extracted  using  impedance-based  techniques  [4] . 

A(s)  s2X(s) 

W)  =  ^W 


where  A(s)  is  the  acceleration  of  the  mass 


and  F(s)  is  the  input  force  to  the  mass. 


From  the  generalized 
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spring-mass-damper  model,  we  have 


F(s)  =  mts2X(s )  +  H(s)X(s ) 

where  the  second  term  represents  the  combined  damping  and  restoring  forces,  i.e. 

■  +  Fd( 

For  a  classical  mechanical  oscillator,  this  term  would  simply  be 

H(s)X(s)  =  ( k  +  sb)X(s ) 


F(s')X(s')  —  FSpring  -f-  r damp 


(2) 

(3) 

(4) 


If  we  wish  to  consider  a  general  form  in  which  the  combined  forces  are  defined  by  a  fractional  derivative,  we 
have 

H(s)X(s)  =  (k  +  sa/3 )  X(s)  (5) 

where  k  is  a  pure  capacitive  term,  (3  is  a  combined  capacitive  and  resistive  term,  and  a  is  the  order  of  the 
fractional  differentiation. 

Returning  to  Equation  2  and  isolating  this  term  yields 


F(s)  H(s) 

—  Ttlt  H - y — 


s2X(s) 


H(s)  =  s 2 


F(s)  1 

2 

'mwv1 

Ux(s)  mt\ 

=  8 

(6) 

(7) 


Where  ^  is  the  inverse  of  the  experimentally  determined  transfer  function,  and  mt  is  the  sum  of  the 
masses  of  the  titanium  cube,  foam  pads,  and  accelerometer,  as  measured  using  a  digital  scale.  Thus,  the 
experimental  data  that  represents  the  bracketed  term  on  the  right  hand  side  of  Equation  7  is  multiplied  by 
s 2  to  yield  the  data  that  represents  H(s),  i.e., 


H(s)  =  -u 


A(s) 

F(s) 


-1 


mt 


(8) 


Note  that  H(s)  in  Equation  8  has  units  of  kg/s2.  Solving  Equation  5  for  H(s )  and  equating  with  Equation 
8  begins  to  illustrate  the  procedure  for  characterizing  the  material  behavior 


H(s)  =  (k  +  sa/3)  =  -LO2 


Ms) 

F(s) 


-1 


mt 


Replacing  the  lumped-parameter  constants  k  and  f3  with  the  following  bulk  material  moduli, 

.-2Ar  n-2Ar 


(9) 


(10) 


where  A  and  S  are  the  foam  pad  cross-sectional  area  and  thickness,  respectively,  results  in  an  equation  that 
relates  the  model  to  experimental  data 


F(s) 


-  mt 


(11) 


L 


J 


model  experimental  data 

This  demonstrates  that  the  experimental  data  can  be  represented  by  a  three-parameter  model  of  the  material 
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behavior:  Go,  G i,  and  a  ,  as  formulated  by  Bagley,  et  al  [1].  Setting  s  =  iuj  on  the  left  side  of  Equation  11 
and  applying  Euler’s  identity  yields  the  following 

(O/'K  Ol'K  \ 

cos  — — j-  i  sin  — J  (12) 

From  Equation  12  ,  the  real  and  imaginary  parts  of  Equation  11  can  be  expressed  as 

Re  (data)  =  Go  +  ^Gi  cos  cja  (13) 

Im  (data)  =  (Gi  sin  uoa  (14) 

A  plot  of  the  real  and  imaginary  parts  of  the  experimentally  determined  transfer  function  for  a  particular 
material  is  shown  in  Figure  7  Where  Re  (data)  and  Im  (data)  are  the  real  and  imaginary  parts  of  the  right  side 


600 


frequency  (Hz) 


Figure  1:  Real  and  Imaginary  parts  of  transfer  function 


of  Equation  11,  as  determined  from  the  experimental  data.  Taking  the  logarithm  of  both  sides  of  Equation 
14  indicates  a  linear  trend  -  where  logo;  is  the  independent  variable  -  from  which  the  order  of  fractional 
differentiation,  <a,  can  be  estimated 


log  [Im  (data)]  =  a  logcj  +  log  ^Gi  sin  —  J  (15) 

Thus,  a  can  be  obtained  from  the  following  using  the  resulting  slope  and  y-intercept  from  a  linear  regression, 
m  and  5,  respectively 

(CX7T  \ 

G i  sin  —)  (16) 

Using  a  base  ten  logarithm,  the  equation  for  Gi  becomes 


Gi 


106 

sin  ^ 


(17) 


Having  obtained  two  of  the  three  parameters  of  the  model,  Go  remains  to  be  determined.  Substituting  the 
expressions  for  a  and  G\  into  equation  13  and  rearranging  yields 


Go  =  10Re 


(18) 


Note  that,  according  to  the  above  equation,  Go  is  a  function  of  the  frequency,  cj,  and  thus  would  not  appear 
to  fit  into  the  scheme  of  a  linear  model.  However,  a  plot  of  Go  vs  a;,  for  a  typical  viscoelastic  foam  shows 
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that  there  is  minimal  variation,  i.e. ,  less  than  three  percent  error,  over  the  frequency  range  of  interest,  as 
indicated  in  Figures  8  and  9  .  A  plot  of  the  curve  fit  for  the  experimentally  determined  transfer  function 


Figure  8:  comparison  of  Gq(lj)  and  constant  Go 


for  a  particular  candidate  material  is  shown  in  Figure  10  The  curve  fit  parameters  are  the  experimentally 
determined  material  properties,  as  given  in  Table  2. 


padding  material 

color 

a 

Go  (N/m2) 

Gi  (N-sa/m2) 

Ethylene  Vinyl  Acetate 

white 

0.56 

3.00  x  105 

3.19  x  103 

Acrylonitrile  Butadiene  Rubber  /  PVC 

black 

0.60 

2.39  x  105 

3.05  x  103 

Table  2:  Experimentally  Determined  Properties  of  Padding  Materials 
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Figure  10:  Curve  fit  for  experimentally  determined  transfer  function 
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Figure  11:  Computational  approach  for  mixed  state-space  equations 


5  Computational  and  Simulation  Techniques 

The  equations  of  state  extracted  from  the  bond  graph  model  constitute  a  mixed  state-space  formulation  of 
the  system  of  differential  equations,  whose  general  form  is  shown  in  Equation  19. 


x  =  A 


[  XI  i 

[foil 

<•••>+  B  u(t)  +  $  D“ 

... 

lXF  J 

LUfJJ 

(19) 


The  solution  requires  an  approach  where  the  fractional  derivative  of  a  state  variable  is  explicitly  computed  [3] . 
Figure  11  illustrates  this  approach,  where  the  differential  equations  are  solved  using  an  existing  variable-step 
fourth-order  Runge-Kutta  algorithm  which  is  linked  to  a  subroutine  that  computes  the  fractional  derivatives 
of  the  necessary  state  variables  only ,  at  each  time  step,  t*.  The  simulations  are  run  using  MATLAB™with 
a  fixed  time-step  size  of  10-5  sec. 


6  Simulation  Results 

To  assess  the  validity  of  the  state  equations  that  model  the  dynamics  of  the  helmet  and  padding  materials, 
simulations  were  run  to  determine  whether  the  simulation  results  would  match  the  experimental  results 
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Elhylene  Vinyl  Acetate  (while  foam) 


Figure  12:  Validation  of  state  equations  via  comparison  with  experimental  data 
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Figure  13:  Summary  of  all  simulation  runs  and  dynamic  behavior  metrics  for  various  blast  intensities  and  padding 

materials 


obtained  as  described  in  Section  3,  i.e. ,  the  input  force  was  set  to  the  that  of  the  impact  hammer  for  the 
material  property  experiments,  instead  of  the  forces  obtained  from  the  CFD  models  of  the  explosions.  It 
was  demonstrated  that  there  was  good  agreement  between  the  simulated  and  experimental  results,  as  shown 
by  the  overlay  plot  in  Figure  ??.  Simulations  were  run  for  both  materials  subjected  to  a  range  of  blast 
intensities,  as  gather  from  the  CFD  analysis  study  cited  in  Section  2.  The  metrics  used  to  identify  material 
behavior  were: 

1.  The  peak  force  transmitted  to  the  victim  by  the  padding  material,  expressed  in  newtons  and  lb-force 

2.  The  attenuation  factor,  i.e.,  the  ratio  of  the  peak  transmitted  force  and  the  peak  blast  force,  expressed 
as  a  percentage 

3.  The  maximum  compression  of  the  padding  material,  expressed  in  mm 

4.  Whether  or  not  the  material  became  saturated,  i.e.,  the  material  reached  its  maximum  possible 
compression 

Figure  12  summarizes  the  simulation  results  for  various  blast  intensities  and  padding  materials.  Example 
simulation  results  for  cases  of  material  non-saturation  and  saturation  are  shown  in  Figures  13  and  14, 
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respectively.  These  plots  and  the  summary  of  results  demonstrate  how  the  combined  experimental  and 
simulation-based  approach  for  assessing  the  efficacy  of  materials  in  blast  applications  for  actual  padding 
materials  that  are  to  be  considered  for  helmet  and  other  PPE  applications. 
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